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STABLE PLASMA CONFINEMENT BY MULTIPOLE FIELDS* 
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Multipole magnetic fields increase with distance 
from and are convex toward the center of sym- 
metry, providing a central equilibrium position 
and correct field curvature for stability of con- 
fined plasma. Braginskii and Kadomtsev' and 
Grad? have treated stability only for this convex 
portion of field, and Tuck® pointed out the utility 
of joining lines of force exterior to the multipole 
conductors to eliminate plasma escape along lines 
of force. 

We describe a multipole configuration for which 
it has been possible to derive the conditions for 
plasma stability over the completely contained 
tubes of flux. These conditions can conveniently 
be met. The structure has only r and @ fields; 
thus eliminating the z field of Stellarators and 
pinches, and critical surfaces which break up 
into multiple pressure-reducing flux bundles un- 
der the influence of field errors.* A large volume 
of high B(=2p/B?) plasma is contained with no 
plasma surfaces in the sense of the cusp model.'~ 
The configuration is stable due to the externally 
supplied fields alone, no matter how little plasma 
is present. There is no requirement for persist- 
ence of special plasma current distribution as 
employed in pinches. 

Figure 1 shows the cross section of a multipole 
field configuration. The z direction can be linear 
or along the circumference of a torus. Toroidal 
effects were found sufficiently small to be omitted 
in the calculations. Current flows in one direction 
through the rods and in the opposite direction on 
the wall. Briefly, the stability of the exterior of 
the multipole results from the flux tube connection 


between Region I where the field is convex toward 
the plasma, thus providing stability, and Region 
II where the field outside the plasma is concave 
toward the plasma contributing instability to the 
flux tube. The plasma on flux lines which circle 
a rod is in a completely stable situation. Condi- 
tions for net stability follow from the considera- 
tions sketched below which will be fully published 
elsewhere. 

The energy principle’ is used with the flux func- 
tion, y, and orthogonal coordinate, x, chosen so 
that 


B,= -8y/8r, 


" (1/r)(8y/26), 
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Cross section of multipole configuration. 


FIG, 1. 
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and so that a line element, ds, is 
(ds)? = (1/B?)(d y)? +J*B?(dx)? + (dz). 


We minimize the energy, 5W, with respect to 
displacements in the x andz directions, elimin- 
ating terms containing these components. The 
sufficient condition for stability that results is 
that the function, 


S=§dy{B “2g *(ax, /8 9)? 
-X,>JB(ap/8y)[0(p + 3B)/8y]}, 


be positive for all flux lines, where Xp is defined 
by 


— = 2) (1/B)X, coskz, 
¥ p=0 ’ 


and by is the ¥ component of the displacement. 
Consider first a uniform X, along the arcs in Fig. 
2. Then the stability condition becomes 


(B,?/B,”) < o,/ 2: 


This condition is easy to satisfy since B, may be 
an order of magnitude less than B,. 

A worse displacement is one which does not 
occur in Region I where stability is provided, but 
rather the displacement is only in Region II. Sta- 
bility would then have to come from bending the 
lines of force. The stability condition for this 
type of displacement is 


(dp /dr,) <(B,?/2r,)(1/ 0,4), 


REGION II 





REGION I 


FIG, 2. Geometry of flux line for stability calculation. 
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or approximately 
B,< (d/r,)(0/ 4). 


This condition can also be satisfied with conven- 
ient proportions of rod and rod cavity. 
A general result for plasma containment time 
is 
o 2 
T~V nax / (MP o)» 
where 7 is plasma resistivity in emu, or 


-10-77 92 2 
T=10 d*(B,/B,) seconds. 


This time is generally long compared to the time 
for other processes to affect the plasma. 

Even at low density and in the collisionless case, 
pressure is scalar and constant along a field line. 
A spherically symmetric velocity distribution in 
the central low field is maintained in the high- 
field region where the orbit magnetic moment is 
an invariant. In case the ion gyroradius is less 
than d, electrostatic ion reflection due to electron 
space charge occurs in the high-field region and 
in this case the action is invariant. This insures 
constancy of pressure on the central field line. 

A free rod can have a stable equilibrium posi- 
tion, but with supports for the toroidal case, plas- 
ma lifetimes useful for some purposes are possi- 
ble. 

Linear arrays allow rods to be connected out- 
side the plasma. Termination in a point cusp re- 
quires the cusp flux to join the rod flux and to 
wind its way along the multipole, producing a B,. 
There is a confluence of cusp flux, as it starts 
around the rod, where a singular line is generated 

Injection of plasma through the convex surface 
of Region I allows entry by the instability caused 
by reversal of the pressure gradient.® Once with- 
in, the plasma is confronted with a stable con- 
figuration which confines it. 
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EXPERIMENTAL EVIDENCE FOR QUANTIZED FLUX IN SUPERCONDUCTING CYLINDERS* 


Bascom S. Deaver, Jr.,and William M. Fairbank 
Department of Physics, Stanford University, Stanford, California 
(Received June 16, 1961) 


We have observed experimentally quantized 
values of magnetic flux trapped in hollow super- 
conducting cylinders. That such an effect might 
occur was originally suggested by London’ and 
Onsager,” the predicted unit being hc/e. The 
quantized unit we find experimentally is not hc/e, 
but hc/2e within experimental error.*® 

Although the unit of quantized flux is small 
(hc /2e = 2.0710 gauss cm’), it can be produced 
by a magnetic field easily measured and con- 
trolled in the laboratory if the area to which it 
is confined is sufficiently small. For our sam- 
ples, one flux unit corresponds to a magnetic 
field of the order of 0.1 gauss. Measurements 
were made on two hollow tin cylinders. Cylinder 
No. 1 was 0.8 cm long, 2.33 x10-* cm outside 
diameter and 1.33 x107° cm inside diameter. 
Cylinder No. 2 was 0.9 cm long, 1.64107 cm 
outside diameter and 1.35x10~* cm inside di- 
ameter. These were fabricated by electroplat- 
ing tin on a one-centimeter portion of a No. 56 
copper wire. The sample, plus protruding wire, 
was jacketed for protection and strength with 
electroplated copper to an approximate outside 
diameter of 8x10-° cm. 

A field-free region (H=0+ 0.001 gauss) is 
prepared using three orthogonal Helmholtz coils. 
The tin cylinder is placed in this region and 
cooled through the superconducting transition 
in the presence of a known applied axial magnetic 
field. The net flux in the cylinder is measured 
both with the field on and after the field is turned 
off. The measurement is made by moving the tin 
cylinder up and down one hundred times per sec- 
ond with an amplitude of one millimeter and ob- 
serving the electrical pickup in two small coils, 
each of ten thousand turns, surrounding the ends 
of the cylinder. The instrument is similar in 
concept to that described by Foner.* The induced 
emf measures the difference in the flux contained 
within the area of the cylinder and that which 
would have been in the same area if the cylinder 
were absent (or in the normal state). The system 
is calibrated by cooling the sample below the 
superconducting transition in zero field and ob- 
serving the signal from the completely diamag - 
netic cylinder when a known magnetic field is 
applied. From the value of the field and the meas- 
ured cross-section area of the cylinder, the ab- 


solute value of the flux for a given signal is cal- 
culated. 

The diameter of each cylinder was measured 
with a microscope equipped with a micrometer 
eyepiece. X-ray photographs verified the di- 
mensions of the tin cylinder after the application 
of the copper jacket. For the purpose of cal- 
culating flux, the measured radii of the cylinders 
are reduced by 0.6 micron due to an expected 
loss of superconducting properties on the sur- 
face of the tin in contact with the copper.® That 
this correction was approximately valid is in- 
dicated by a 0.2°K decrease in the value of the 
transition temperature for the sample No. 2 
whose cylindrical walls were 1.5 microns thick, 
leaving, we believe, only 0.3 micron of super- 
conducting material in the center after allowance 
for the effect of the center copper wire and the 
outside copper jacket. This is in agreement with 
the experimental results on electroplated tin.® 

With this adjustment, the area used for the 
diamagnetic calibration of cylinder No. 1 is 
3.84 x10~* cm’, and the area of the hole is 
1.65x10~* cm?. For cylinder No. 2 the corres- 
ponding quantities are 1.81 x10-® cm? and 1.70 
x107* cm?. 

Data on sample No. 1 are shown in Fig. 1, and 
on sample No. 2 in Fig. 2. The diagonal line 
through the origin represents the calibration. 

It is the signal corresponding to zero flux in the 
cylinder and hole in the presence of the applied 
field as described above. The experimental 
points shown on the graph represent two types of 
data for each value of the applied field. The 
points on the bottom half of each graph represent 
the signal in the presence of the applied field 
after cooling through the transition in that field. 
The points in the upper half represent the trapped 
flux after the field is subsequently reduced to 
zero. For each point in the lower curve there 

is a corresponding point in the upper curve. The 
solid lines represent calculated integral values 
of hc /2e. 

It can be seen that certain features of the data 
are common to both samples. (1) Below a cer- 
tain value of applied field, the total cross section 
of the cylinder acts as a perfect diamagnet, ex- 
cluding all the flux, and no flux is trapped when 
the applied field is turned off. (We believe this 
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FIG. 1. (Upper) Trapped flux 
° 9° in cylinder No. 1 as a function 
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of magnetic field in which the 
cylinder was cooled below the 
superconducting transition tem- 
perature. The open circles are 
individual data points. The solid 
circles represent the average 
value of all data points at a par- 
ticular value of applied field in- 
cluding all the points plotted 
and additional data which could 
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two hundred data points are rep- 
resented. The lines are drawn 
at multiples of hc/2e. (Lower) 
Net flux in cylinder No. 1 be- 
fore turning off the applied field 
in which it was cooled as a func- 
tion of the applied field. Open 
and solid circles have the same 
significance as above. The low- 
er line is the diamagnetic cali- 
bration to which all runs have 
been normalized. The other lines 
are translated vertically by suc- 
cessive steps of hc/2e. 








provides a way of obtaining a truly zero-magnetic - 
field region.) (2) When the applied field exceeds 

a certain value, flux is trapped both with the field 
on, and after the applied field is turned off. The 
amount of this trapped flux within the experi- 
mental accuracy of the data is hce/2e. 

The amount of trapped flux remains constant 
as the applied field is increased until a value 
approximately three times that for the initial 
trapping is reached, at which point the trapped 
flux increases to about twice the original amount. 
There appears to be evidence for additional 
changes at five and seven times the field for the 
first trapping. 

Fluctuations in the data are caused by vari- 
ations in the zero of the magnetic field, changes 
in the gain, vibration amplitude, drift, and ran- 
dom noise in the detection system. The approxi- 
mately two hundred data points for sample No. 1 
were taken over a three-week period during 
which the drift and noise were gradually im- 
proved. The fluctuations of the data around the 
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values 0 andhc/2e represent, we believe, ex- 
pected scatter from drift and noise. This scat- 
ter has been greatly improved for sample No. 2. 
For both samples the data are consistent with 
values 0 and hc/2e for the trapped flux as de- 
scribed above. 

Near the transition to the second and third 
steps the fluctuations in the data are greater, 
and in addition points lie between the steps. 
Some increased scatter is expected since the 
absolute fluctuations due to changes in gain and 
vibration amplitude are proportional to the size 
of the signal. The points between the steps do 
not necessarily indicate trapping of nonintegral 
values of flux. Since the observed signal is the 
sum of the emf’s from coils at the two ends of 
the sample, a flux line passing out of the cylin- 
der at some point other than the end may pro- 
duce different signals in the two coils. The two 
ends of the cylinder are not quite identical; so 
near the transition region it is probable that the 
two ends might trap a different number of units 
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of flux, the extra unit being shoved out the side 
of the cylinder. This is especially probable for 
sample No. 1 since the x-ray photograph showed 
a break in the tin coating near the middle of the 
cylinder. Also, it is known that flux can create 

a normal region in a superconductor by shrink- 
ing in size until the critical field is exceeded. 

In this experiment we were unable to measure 
independently the signals from the two coils. 
However, in future experiments this will be done 
to remove this ambiguity. It is interesting to note 
that no intermediate points are found outside the 
expected scatter of the data near the first step. 
One point for which no flux was trapped was 

found near the center of the first step with sample 
No. 1. 

In conclusion, we find: 

1. The flux trapped in a superconducting cylin- 
der both in the presence and absence of an applied 
magnetic field is not continuous but exhibits a 
step behavior, the first step occurring for 
@=hc/2e, within experimental error in the data. 
Considering all sources of error,, we estimate 





that the value of the trapped flux at the first 
step is hc/2e+20%. If the correction to the 
size of the cylinder due to the presence of the 
copper should prove invalid, an additional 11% 
error could arise for the large cylinder and 17% 
for the small cylinder. 

2. The data seem to indicate additional steps 
athc/e, 3hc/2e, and 2hc/e. The points appear- 
ing between these levels will be investigated 
further. 

3. The ratio of the fields at which the steps 
occur are approximately 1, 3, 5, and 7. In the 
first cylinder (for which the effective cross-sec- 
tional area of the cylinder is 2.33 times the area 
of the hole), the first jump occurs when the flux 
passing through the total effective cross section 
of the cylinder in the normal state is approxi- 
mately hc/2e. 

For cylinder No. 2 (in which the effective cross- 
sectional area of the cylinder is 1.1 times the 
area of the hole), the first jump occurs when the 
flux passing through the total effective cross sec- 
tion of the cylinder in the normal state is approxi- 
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mately 0.6hc/2e. In a following Letter, Byers 
and Yang® conclude that in a thin ring the first 
jump should occur at 0.5hc/2e. 

4. Since the time constant of our measuring 
circuit is 25 seconds, this experiment gives 
only a large upper limit for the time involved 
in reaching these quantized flux values. Mer- 
cereau and Vant-Hull’ have reported a negative 
experiment designed to observe quantized flux 
in a l-mm ring cooled 6000 times per second 
through the superconducting transition in a small 
magnetic field. It is possible that the difference 
in their results and the results of our experiment 
are due to a minimum time necessary to estab- 
lish equilibrium. We are planning to investigate 
this relaxation time. 

We have had the pleasure of discussing the 
results of this experiment with N. Byers, C. N. 
Yang, and L. Onsager, whose interpretation of 
these results appear in the following Letters. °® 
One of us (WMF) also wishes to acknowledge his 
indebtedness to F. London and M. J. Buckingham 
who greatly influenced his concept of the super- 
fluid state. We also wish to thank F. Bloch, L. L 
Schiff, and J. D. Bjorken for many stimulating 


discussions of the experiment. We wish to ac- 
knowledge the invaluable assistance of M. B. 
Goodwin. 
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’ THEORETICAL CONSIDERATIONS CONCERNING QUANTIZED MAGNETIC FLUX 
IN SUPERCONDUCTING CYLINDERS* 
N. Byers and C. N. Yangt 


Institute of Theoretical Physics, Department of Physics, Stanford University, Stanford, California 
(Received June 16, 1961) 


In a recent experiment,! the magnetic flux 
through a superconducting ring has been found 
to be quantized in units of ch/2e. Quantization 
in twice this unit has been briefly discussed by 
London? and by Onsager.* Onsager‘ has also 
considered the possibility of quantization in 
units ch/2e due to pairs of electrons forming 
quasi-bosons. 

The previous discussions® leave unresolved 
the question whether quantization of the flux is 
a new physical principle or not. Furthermore, 
sometimes the discussions seem? to be based on 
the assumption that the wave function of the su- 
perconductor in the presence of the flux is pro- 
portional to that in its absence, an assumption 
which is not correct. We shall show in this Let- 
ter that (i) no new physical principle is involved 
in the requirement of the quantization of magnetic 
flux through a superconducting ring, (ii) the 
Meissner effect is closely related to the require- 
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ment that the flux through any area with a bound- 
ary lying entirely in superconductors is quan- 
tized, and (iii) the quantization of flux is an 
indication of the pairing of the electrons in the 
superconductor. 

Macroscopic discussion. Consider a multiply 
connected superconducting body P with a tunnel 
O (Fig. 1). We shall only discuss macroscopic 








FIG, 1. Multiply connected superconductor. 
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dimensions much larger than the penetration 
depth. The Meissner effect then states that in- 
side the superconductor P the magnetic field is 
zero, and the current is zero. Surface currents, 
however, do exist and persist on the surfaces S, 
and S,. The surface currents and the external 
sources of magnetic fields together produce no 
magnetic flux in the interior P of the supercon- 
ductor. They in general, however, produce a 
net magnetic flux through O, to be denoted by ®. 

The energy eigenfunction y of the electrons in 
the superconductor satisfies 


1 > Cr 
2a Fog LV +A Py + Vp =Ey, (1) 


J 


where A is the vector potential due to the sur- 
face currents and external magnetic sources. 
Inside P, 


vxA=0. 


Hence A = Vx, where x is not single valued in P 
but increases by 


Ay =fA-di = JJH-do = 4, (2) 


whenever one goes around the tunnel O once. De- 
fining 


‘— ; é = 
v) =vexpLigxlt ), (3) 
J 
we see that (1) reduces to 


1 o , i , 
Lam (ie vy +Vy'=Ey’. (4) 


The vector potential A is eliminated from this 
equation. However, the boundary condition for 
vy’ is that when all electron coordinates are 
fixed, except for one, f;, and fF, is brought 
around O once, y~’ changes by a constant factor 
yaya (Me, (5) 


[To prove (5) we use (3) and (2) and the fact that 
~ is single valued. | 

The eigenvalues E are determined by the dif- 
ferential equation (4) and the boundary condition 
(5). It is thus obvious that we have: 

Theorem 1. The energy levels are periodic 
in the magnetic flux @ with a period ch/e. 

If the surfaces S, and S, are concentric cylin- 
ders on which ~=0, and V is put equal to zero, 
the energy levels E can be explicitly solved for, 
illustrating this theorem. One notices that y’ is 
not simply proportional to y(@=0), as is some- 
times? assumed ‘in, the literature. 
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FIG. 2. Periodic variations in N~'1nQ as a function 
of trapped flux ®, 


If V is a real function of Yj» by taking the com- 
plex conjugate of (4) and (5), we have: 

Theorem 2. The energy levels are even func- 
tions of #. 

It is clear that theorems 1 and 2 remain valid 
if we introduce the lattice coordinates of the 
metal and if we introduce the spin of the elec- 
trons. The operation of complex conjugation in 
the proof of theorem 2 has then to be replaced 
by the time reversal operation and the proof 
depends on the time reversal invariance of the 
interactions. 

From these theorems it follows trivially that: 

Theorem 3. The partition function Q of the 
system is an even periodic function of @ with 
period ch/e. 

At @=(ch/2e xinteger, this theorem shows that 


alng /ae@ =0, (6) 


and that lnQ@ has the general form shown in Fig. 2, 


Now the body current in the superconductor 
around O is 


1=kTc8 1nQ/2®. (7) 


(c =velocity of light.) In the differentiation we 
keep the temperature T constant. 

The Meissner effect requires that J/=0. Thus 
the equilibrium states are given by the maxima 
and minima on the curve in Fig. 2. We shall 
now give an argument to show that the maxima, 
not the minima, are the equilibrium states real- 
ized. A point D in Fig. 2 is not an equilibrium 
state so the calculation of the partition function 
at that point is strictly speaking meaningless. 
But the slope of the curve at that point indicates 
that if a flux 4, is made to pass through O, a 
body current would be induced, the sense of 
the current being negative according to (7). 

The additional flux due to this body current 
causes the flux through O to decrease. The 
equilibr'um state reached would therefore be E 
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which is a maximum of the curve. We state this 
as: 

Theorem 4. The superconducting state is 
given by the maxima of In@ as a function of 4. 

If the external flux does not assume a value for 
which In@ is a maximum, surface currents will 
flow on S, and S, to make up a total flux @ for 
which InQ is a maximum. 

The experiment of reference 1 and theorem 4 
together prove that lIn@ has maxima at integral 
values of @/(ch/2e). Whether a microscopic 
theory yields these maxima will be discussed 
in the next section. If it does, then theorem 4 
shows the following: The flux through any sur- 
face whose boundary loop lies entirely in super- 
conductors is quantized in units ch/2e. The re- 
quirement of this quantization in turn clearly im- 
plies that the flux through any small area ina 
superconductor is zero; hence it implies the 
Meissner effect. 

In a loose sense the above argument can be 
used to “derive” the Meissner effect itself: If 
the magnetic flux in a superconductor is not zero, 
body currents will flow around all loops through 
which the flux is not quantized. The system can- 
not reach a steady state until all magnetic flux 
is expelled from the interior of the superconduc- 
tor. 

Microscopic considerations. We now want to 
see whether a microscopic calculation does or 
does not lead to maxima of InQ at @/(ch/2e) 
= integer. 

To investigate this point we first take a collec- 
tion of noninteracting spinless electrons between 
two concentric cylinders S, and S,. The electrons 
at a point at a distance ry from the axis have mo- 
menta p,, pg, and pz in the radial, azimuthal, 
and z directions. Clearly 





Dot =nh. (n=integer) 


The energy of the electron is 


2 2 

ben Seay} 0 
The partition function Q can be computed from 
such an energy spectrum. The resultant N~'InQ 
[to the order N°] does not depend on #. Thus, 
according to theorem 4, for a collection of non- 
interacting electrons, the flux in the tunnel does 
not have to be quantized. 

It is not difficult to understand why for such a 
model N~1nQ@ does not depend on ®. To see this 
we suppress the p, and p, degrees of freedom 
and take the temperature T7=0. The one-dimen- 
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sional Fermi sea problem, 


hi? e ¥ 
2mr Dfe th ‘) , 
gives an average energy per particle of 


E ee 
anal ——_}| 
N constant + Imrcne® ’ IF | 2) (9) 


if N =the number of particles is odd. But if N is 
even, then 


E_ h? ® ,V 1 eo 
x7 constant + =——a| 4) ._ = << - iate (10) 
hk? ko ed 
= constant tee} . 0> 7 > -2. 
(10’) 


Thus, depending on the evenness or oddness of N, 
the energy has a minimum at e@/ch =4 or 0 (mod- 
ulo 1). The three-dimensional problem at T=0 
is decomposable into many one-dimensional prob- 
lems with varying values of N. Thus the above- 
discussed fluctuation leads to a cancellation for 
the three-dimensional problem, resulting in an 
N~lnQ versus @ curve that is flat. (An N“lnQ 
curve that is flat applies to the case of a metal 

in its normal rather than superconducting state.) 
A similar cancellation obtains for T +0. 

In the neighborhood of @=+0 the states with n 
>0 have energies that increase with $, and those 
with n<0 have energies that decrease with ©. 

The average energy for the two states n and -n, 
however, increases with @ like 


constant + (f?/2m7r)(e@/ch)?. (11) 


If there is a “pair correlation” of the kind pro- 
posed by Bardeen, Cooper, and Schrieffer® for 
the superconductor so that states m and -n (ora 
pair of time-reversed states) are either both oc- 
cupied or both unoccupied, (11) becomes the cor- 
rect energy per particle for small @. (In sucha 
case the fluctuation and cancellation phenomena 
disappear.) This is represented in Fig. 3 by the 
parabolas at 2e@/ch = -2,0, 2, etc. 


N“inQ 


=2 -| oO | _ 2eo 
~ ch 





AA NAINKT NY 


FIG, 3. A curvecf N~1nQ versus 2e@/ch, showing 
parabolic behavior near maxima at 2e¢/ch = integer. 
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At 2e@/ch=1, pairing between 


e 


aoe =" + aad and -3, 3 and -3, etc., 


n+ 
occurs and the energy per particle remains the 
same as for the case #=0 (to order N°). In the 
neighborhood of 2e@/ch =1, the additional energy 
for each of these pairs is again twice 


(h? /2mr*)[(e&/ch) -2P, 


which give rise to the parabolas at 2eé/ch =+1 
in Fig. 3. In the absence of a detailed theory, 
we draw a smooth curve in Fig. 3 to extrapolate 
between the maxima. 

Thus the Bardeen, Cooper, Schrieffer pairs 
for the superconducting state give rise to the 
curves such as those depicted in Fig. 3, where 
the parabolas are repeated at periods A(2e@/ch) 
=1, and the central parabola is given by 


hi? ea) 
N~inQ= - fant —(eh) +constant, (12) 
kT2m( laws ch 


where f =fraction of electrons that are paired. 
It is interesting to estimate the magnitude of 
the body current at, say, 


0<2e6/ch<3. (13) 
It is, by (7), 
[= -Nfc(e*/mc*)®/(An*{r*) |). 


The flux induced by this current is, for a thin 
ring superconductor, 


} ait fx(number of electrons in a 


length e?/mc’ of the ring)®. 


For the experiment of reference 1, -®jnduced/® 
>» 1 if f is not too small, showing that the maxima 
in Fig. 3 are very pronounced. 

From Fig. 3 and the argument preceding theo- 
rem 4, we conclude that the trapped flux @ and 


the original flux #, are related in the following 
way: 





,/(ch/2e) } 
4-4 0 
a>% 1 
i> 2 
etc. 


It is interesting to notice that the existence of 
the variation of the energy levels of the electrons 
in P with the flux @, even when there is no mag- 
netic field in P, is based on the same principle 
as the experiment proposed by Aharonov and 
Bohm.*” 

We wish to thank W. M. Fairbank and B. Deaver 
for informing us of the progress of their beautiful 
experiment and for many discussions. We also 
wish to thank F. Bloch for stimulating discussions. 
One of us (CNY) takes this opportunity to thank the 
members of the Physics Department of Stanford 
University for the hospitality extended him during 
his visit. 





*Work supported in part by the U. S. Air Force 
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MAGNETIC FLUX THROUGH A SUPERCONDUCTING RING 






Lars Onsager* 


University of California, La Jolla, California 
(Received June 15, 1961) 


London’ recognized that the magnetic flux 
embraced by a superconducting ring ought to be 
quantized. He argued as follows: The current 
density is proportional to the average of 


p -(Q/c)A, 


where A denotes the vector potential and Q the 
charge of the current-carrying particle. In the 
interior of the superconductor the current den- 
sity vanishes, so that the condition 


$ P+ ax= -nh 
for a single-valued wave function implies 
=F K- dx= -nhc/Q. 


Substituting Q = -e for the charge of the electron, 
he arrived at the conclusion 
@=nhc/e. 


It is possible to recast London’s discussion in 
a form which is completely Lorentz and gauge 
invariant; the details need not concern us here. 

London’s result inspired the suggestion? that 
the quantization of flux might be an intrinsic 
property of the electromagnetic field. 

Not much later, Schafroth*® pointed out that 
electron pairs held together by attractive inter- 
actions would obey Bose statistics and be capable 
of superfluid (Einstein) condensation. A likely 
source of the requisite attractive interaction— 
by way of the phonon field—had been suggested 
by Frdhlich* and by Bardeen.® These ideas form 
the basis of more detailed theories, which ex- 
plain the various observed properties of super- 
conductors so well that they have been generally 
accepted. 

Deaver and Fairbank® have found that the flux 
embraced by a superconducting annulus varies 
in steps of half the size proposed by London. 
This is readily explained by Schafroth’s theory, 
which requires 

Q= 2e, 
@=nhc/2e. 
The discovery of steps just this size provides 





a convincing direct proof for the pairing of elec- 
trons. 

The notion that the electromagnetic field it- 
self might be subject to a similar condition seems 
untenable now, for singly charged bosons exist 
(deuterons) and a condition imposed on the elec- 
tromagnetic field ought to be equally compatible 
with all charged particles. 

Instead, we arrive at the remarkable result 
that we can measure the magnetic flux (except 
for an additive undetermined multiple of hc/ 2e) 
embraced by a given closed path without examin- 
ing the enclosed field; a superconductor placed 
along the path will respond with a supercurrent 
which compensates the fractional excess of flux. 
Complete Meissner effect in a multiply-connected 
superconductor requires coherent ring closure. 
We may infer that such closure will not take place 
unless the free energy liberated by the matching 
of phases exceeds the kinetic energy of the neces- 
sary supercurrent plus the added magnetic field 
energy. The detailed kinetic mechanism is not 
yet known. 

I am indebted to F. Bloch for a discussion 
which clarified the closure problem, and to 
B. Deaver and W. M. Fairbank for the communi- 
cation of their unpublished results. 





*Permanent address: Sterling Chemistry Laboratory, 
Yale University, New Haven, Connecticut. 
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EXPERIMENTAL PROOF OF MAGNETIC FLUX QUANTIZATION IN A SUPERCONDUCTING RING* 


R. Doll and M. Nabauer 
Kommission fiir Tieftemperaturforschung der Bayerischen Akademie der Wissenschaften, 


Herrsching/Ammersee, Germany 
(Received June 19, 1961) 


From theoretical considerations, based on wave 
mechanics, London’ concluded that the magnetic 
flux frozen in in a twofold-connected superconduct- 
ing body (ring or tube) should not have any arbi- 
trary value, but only such values which are inte- 
ger multiples of a basic unit ¢,, 


$,=hce/e =4.12x10™ gauss cm?, (1) 


That means the magnetic flux should be quantized. 
Bardeen and Schrieffer” also agreed with this con- 
clusion. 
In order to verify a possible flux quantization, 
the mechanical torque exerted by a magnetic 
field H, on a small superconducting lead tube 
with frozen-in magnetic flux has been measured. 
From the following equation, one obtains the 
magnetic field H go necessary to freeze in just 
one flux unit ¢,;: 


Hy =? ox 4/nd’), (2) 


where d is the diameter of the tube. The diameter 
has to be chosen sufficiently small that the fluc- 
tuations of the earth’s magnetic field can be ne- 
glected relative to the magnetic field adequate to 
freeze in one flux unit. This has been achieved 
with a tube of 10.3-micron diameter and 0.6-mm 
length (see Fig. 1), which gives a magnetic field 
of Hy9=0.5 oe, according to Eq. (2). 

The sample consists of a small lead cylinder, 
prepared by evaporating lead on a quartz fiber 
(of about 10-micron diameter and with a length 
of about 1 mm). 





Hy 
FIG, 1. Schematic diagram of the sample with the 


directions of the applied field Hy to be frozen in, and 
the measuring field H,. 


The very small torque, proportional to the fro- 
zen-in flux and to the measuring field H, (normal 
to Hy), can be observed by an already known auto- 
resonance method.** The sample is suspended 
on a thin torsion fiber (normal to Hy and H,) in- 
side a coil. In connection with a mirror for re- 
cording the oscillation amplitude it represents 
a system of damped torsion oscillation. The 
oscillation of the system can be kept at a con- 
stant amplitude by the alternating torque, caused 
by periodically reversing the magnetic field H,. 
of the coil. The oscillating system itself controls 
the time of switching the field by means of a pho- 
toelectric device. The damping of the system 
being known, the constant resonance amplitude 
is a measure for the torque which acts on the 
sample. 

Each experimental value was obtained in the fol- 
lowing manner: 1. The sample was heated above 
the transition temperature; then a defined field to 
be frozen in, Hy, was applied. 2. After recooling 
below the transition temperature, H,, was switched 
off. 3. The resonance amplitude was measured 
as described above. 

The resonance amplitude is proportional to the 
product (measuring field H,.) times (magnetic mo- 
ment of the sample). The latter in turn is propor- 
tional to the frozen-in flux. Figure 2 shows the 
measured resonance amplitude divided by the 
driving field H, as a function of the field Hy. (H,. 
has always been about 10 oe.) 

As Fig. 2 shows, it is impossible to freeze in 
any flux between H, =-0.1 and +0.1 oe. Near +0.1 
oe there occur marked steps. Upon increasing the 
magnetic field H,,, the frozen-in fluxes remain 
nearly constant between 0.1 and 0.3 oe. At 0.3 oe 
another step occurs, again followed by a series of 
constant values. 

This is exactly what is expected of a quantized 
magnetic flux in a twofold-connected superconduct- 
ing body. If an arbitrary flux could be frozen in, 
the relation between magnetic flux and field H 
would be as shown in Fig. 2 by the dashed line. 
This has been obtained by measurements at com- 
paratively high fields (Hy =10 oe), in which case 
the value of one flux unit is already very small 
compared with the entire frozen-in flux. 

With the microscopically measured diameter of 
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FIG, 2. Resonance amplitude divided by measuring 
field H, as a function of the applied field H,. The 
ordinate is proportional to the frozen-in flux. X*-— First 
run; o— second run. 


the lead tube, Eq. (2) predicts for the interval of 
the magnetic field strength corresponding to one 
flux unit a value of H, =0.5 oe. The experimen- 
tally observed interval, however, reaches only 


0.2 oe, that is about 40% of the calculated value. 


So far the reason for this discrepancy is not clear, 


For example, an error of 60% in the determina- 
tion of the lead tube’s diameter would explain the 
difference, but such an error is improbable. 

The experiments are being continued with higher 
fields H. and other superconductors of various di- 
ameters. 

Mercereau and Vant-Hull® also tried to verify 
London’s postulate of the quantization of magnetic 
flux in a superconducting ring. The result of their 
experiments was negative. 


The authors are indebted to Professor W. Meiss- 


ner who made possible and promoted this work. 
The authors would further like to thank Professor 
F. X. Eder for encouragement and helpful discus- 
sions. 
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EXPERIMENTAL DETECTION OF TRANSITION RADIATION*t 


H. Boersch, C. Radeloff, and G. Sauerbrey 


I. Physikalisches Institut of the Technische Universitat Berlin, Berlin, Germany 
(Received February 17, 1961) 


If an electron approaches the boundary between 
vacuum and metal, a changing dipole field due to 
the electron and its image charge is formed 
which, according to the theory of Ginsburg and 
Frank,’ effects the emission of “transition ra- 
diation.” In the experiment electrons bombard- 
ing a metal surface generate a visible radiation 
known as “Lilienfeld radiation.” The results of 
the experiments done before this work are, how- 
ever, contradictory.7* Therefore we undertook 
to exclude by more careful experiments the in- 
fluence of surface contaminations and to com- 
pare the experimental Lilienfeld radiation with 
the theoretical transition radiation by variation 


52 


of different parameters. 

Our investigations were carried out at pres- 
sures of about 10~° mm Hg with massive heat- 
able targets and with condensed films of more 
than 1 yw thickness. The electron beam was pro- 
duced by field emission and pulsed with a fre- 
quency of 100 kc/sec. The energy of the elec- 
trons was in general 2-12 kev. The intensity 
of the Lilienfeld radiation was measured with 
a photomultiplier tube and a phase discrimina- 
tor. The results of our investigations are as 
follows: 

(1) The intensity of the Lilienfeld radiation 
is independent of the temperature of the target. 
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The temperature range was 300°K to 2300°K for 
W targets, up to 1800°K for Ta and Mo, and up 
to 1300°K for Ti. 

The intensity of the Lilienfeld radiation is also 
independent of the gas pressure, as investigated 
for tungsten in the range from 10~° to 10° mm 
Hg. Consequently, if gas films are adsorbed 
on the surfaces, they have no influence. There- 
fore the Lilienfeld radiation is a property of the 
clean metal surface only, like the transition ra- 
diation. 

(2) Like the transition radiation, the Lilienfeld 
radiation of W, Ta, Mo, Ti, and Pt is nearly 
plane polarized (Table I). The electric field 
vector is oscillating in the plane given by the 
normal to the target plane and the direction of 
observation. Its direction is independent of the 
angle of incidence of the electron beam. The 
degree of polarization for Ag and Cs is lower 
than that for the former metals. This decrease 
is probably caused by greater roughness of the 
surfaces. 

(3) The intensities of the transition radiation 
and of the Lilienfeld radiation for W, Ni, Al, and 
Cu increase linearly with both electron energy 
and current density. The investigations ranged 
for electron energy from 6 to 30 kev and for cur- 
rent density from 0.1 to 10 ma/cm”. 

(4) The angular distributions of the Lilienfeld 
radiation for tungsten and of the transition ra- 
diation for a metal with a dielectric constant 
le| =3 are nearly identical from 6 =10° to 6 =60° 
(Fig. 1). 

(5) The spectrum of the Lilienfeld radiation, 


Table I. Degree of polarization P and efficiency na, 
of the Lilienfeld radiation and of the transition radia- 
tion. Angle of observation 9=65° against the normal 
of the target plane. The efficiency 7a, is the ratio 
between the total power of the Lilienfeld radiation in 
the spectral range A =4000 to 5120 A, 4A=1120 A, and 
the power of the electron beam. 











Lilienfeld Transition 

radiation radiation 

(exp.) (theor :) 

Metal P 10 Man P 10 Ay 

Tungsten 0.96 4 1.0 4.3 
Tantalum 0.96 3.5 1.0 3.9 
Molybdenum 0.96 4 1.0 4.3 
Titanium 0.88 2 1.0 3.1 
Platinum 0.92 5 1.0 5.1 
Silver 0.6 8 1.0 5.7 
Cesium 0.5 11 1.0 3.9 





recorded by a quartz spectrograph with a reso- 
lution AA~100 A, forms a continuum (Fig. 2). 
The spectra of the Lilienfeld radiation and of 

the transition radiation for silver are plotted 

in Fig. 3. Both spectra nearly agree in the slow 
increase from red to blue and in the following 
abrupt decrease in the ultraviolet. The theoret- 
ical increase of the transition radiation at 3000 A 
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FIG. 1. Intensity of the investigated Lilienfeld radia- 
tion (solid line) for tungsten and of the calculated tran- 
sition radiation (dashed line) versus the angle of obser- 
vation. «=dielectric constant. 
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FIG, 2. Spectrum of the Lilienfeld radiation for Ta, 
W, Mo, and Ti. Absolute intensity relative to 1-watt 
electron beam power and to a range of wavelength 4A 
=1A. 
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FIG, 3. Spectra of the Lilienfeld radiation (exp.) 
and the transition radiation (theor.) for silver. 


is not reproduced by the experimental spectrum. 
It is probable that in this part of the spectrum 
the sensitivity of our apparatus is not sufficient. 

(6) The efficiency 7 A Of the Lilienfeld radia- 
tion for W, Ta, Mo, Ti, and Pt agrees surpris- 
ingly well with the calculated efficiency of the 
transition radiation (Table I). The efficiency 
for the condensed metals Ag and Cs is only of 
the same order. 

On account of these corresponding character - 
istics the Lilienfeld radiation, produced by met- 
als in the energy range mentioned above, is to 
be interpreted as the transition radiation due to 
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Ginsburg and Frank. We think it impossible to 
interpret the Lilienfeld radiation as the visible 
part of the bremsstrahlung. The Lilienfeld ra- 
diation also cannot be interpreted as “plasma 
radiation” according to the theory of Ferrell,® 
because the plasma radiation only appears on 
foils with two boundary surfaces, whereas the 
investigated Lilienfeld radiation is observed on 
compact metals with only one boundary surface 
within the depth of penetration of the electrons. 

The transition radiation may also contribute 
to the background of the radiation from thin 
foils (two boundary surfaces) which had been 
experimentally studied by Steinmann,® and Brown 
et al.” 





*A detailed publication is intended for the Zeitschrift 


fiir Physik. 
tThis work was supported by the Research Corpora- 


tion, New York. 
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ELECTRONIC EFFECT IN THE ELASTIC CONSTANTS OF GERMANIUM 


L. J. Bruner 


International Business Machines Watson Laboratory, New York, New York 


and 


Robert W. Keyes 
International Business Machines Research Center, Yorktown Heights, New York 
(Received May 25, 1961) 


We have measured the elastic constants of pure 
and heavily doped n-type germanium and discov- 
ered a difference which can be attributed to the 
electronic contribution to the strain energy func- 
tion. The importance of electronic energy bands 
in the theory of the elastic constants of metals 
has been appreciated for a long time. Satisfactory 
quantitative verification of the electronic theories 
of elastic constants has not been possible in metals, 
however, because pertinent parameters of the 
band structures and their dependence on strain are 
not known. Such parameters are known for certain 
semiconductors, however, and, in particular, are 
far more accurately known for germanium than for 
any other solid. Therefore germanium provides 
a unique opportunity for a study of the role of elec- 
tronic energy in the theory of elastic constants. 

We have measured the elastic constants of two 
samples of germanium. One, “pure, ” contains 
less than 10** carriers per cm*. The other, “heav- 
ily doped, ” contains 3.5x10!* arsenic donors per 
cm*’, The heavily doped sample was a }-inch cube 
with exceptionally uniform resistivity, as meas- 
ured by surface probes.’ The measurements 
consisted of the determination of the velocities 
of 10-Mc/sec ultrasonic waves for an appropriate 
set of crystal directions by a conventional pulse- 
echo technique. Measurements were made through 
the temperature range 4.2°K-240°K. The results 
at very low temperatures for the pure specimen 
and the specimen degenerately doped with arsenic 
are given in Table I. The values of C,, as a func- 
tion of temperature are given in Fig. 1. Our values 
of C,, for the pure specimen differ from those of 
McSkimin? by almost 1 % and from those of Fine® 
by 0.5% in the temperature range covered by both 
sets of measurements. 

Referring to Table I, $(C,,+2C,,) and $(C,, -C,,) 
are not significantly different in the two specimens, 
while C,, is 5.5% smaller in the heavily doped 
specimen. Measurements on another specimen 
of similar doping level but not as homogeneous 


confirmed this difference in C,,. The fact that 
only C,, is affected by the addition of the arsenic 


characterizes the difference between the samples 


Table I, The elastic constants of germanium at 
4, 2°K (in units of 10‘! dyne/cm?). 








Specimen Cu (C44 -C 49)/2 (C44 +2C4,)/3 
Pure 6.80 4.06 7.58 
Heavily doped 6.42 4.04 7.66 








as an electronic effect. The terms in the strain 
energy function of the crystal which are propor- 
tional to C,, involve just those components of 
strain which alter the conduction band of germani- 
um by removing the degeneracy of the valleys. On 
the other hand, the terms in the strain energy 
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FIG. 1. The temperature dependence of Cy, in pure 
and in heavily doped n-type germanium. 
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function which are proportional to C,, and C,, con- 
tain only components of strain which do not affect 
the structure of the conduction band of Ge. It can 
be concluded that the change of elastic constants 
results from the addition of electrons to the con- 
duction band. 

The effect on C,, is very large. The addition of 
0.08% impurity changes C,, by 5.5%, or seventy 
times the impurity concentration. The large mag- 
nitude of the change is also in agreement with the 
predictions of the electronic interpretation. We 
have used the usual multivalley-deformation po- 
tential model to construct a theory of the electronic 
contribution to the strain energy function. The 
theory resembles that of Leigh,* who calculated the 
elastic constants of aluminum from a multivalley- 
type model, but is not restricted to the degenerate 
limit assumed by Leigh. According to this model 
the effect arises in the following way: When the 
crystal is strained by a pure shear, the energy of 
some valleys is increased and the energy of other 
is decreased. The mean energy of the valleys is 
constant, so that the electronic energy is not 
changed to first order in the strain. However, in 
the strained crystal the energy can be lowered by 
transferring electrons from the higher energy val- 
leys to the lower energy valleys. This re-estab- 
lishment of statistical equilibrium among the elec- 
trons lowers the energy in terms of second order 
in the strain. Therefore it amounts to a decrease 
in the shear elastic constant. The quantitative 
theory shows that in the low-temperature limit, 


(8C 44/Caa)o= -(4/3)990*9m *E,,2NY8/h7C,,. (1) 


Here m* is the density-of-states mass for one 
valley and the remaining notation is standard. 
For our heavily doped sample (N= 3.5x10'® cm™), 
Eq. (1) gives (6C,,/C,,),= -0.075. The observed 
effect is about three-fourths of this. 

The highest temperature of our measurements 
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is about one-third of the electronic degeneracy 
temperature in the doped crystal. The theory 
predicts that (5C,,/C,,) should increase slightly 
(about 10%) in this range. Figure 1 shows that 
in fact (6C,,/C,,) decreases by about 20%. 

Thus the change in elastic constants produced 
by the solution of arsenic in germanium agrees 
with an electronic interpretation in two respects: 
(1) Only C,, is changed; (2) the magnitude of the 
change is very large, and is given semiquantita- 
tively by the theory. The discrepancies with the 
theory, namely, the difference between the cal- 
culated and observed magnitude of 5C,,/C,, and 
the difference between its observed and predicted 
temperature dependence, probably indicate an 
inadequacy of the theoretical model. However, 
it is difficult to completely rule out any effects 
of inhomogeneities, which are so common in 
heavily doped germanium. Note that the period 
of our ultrasonic wave (10~ sec) is very long 
compared to the scattering times of electrons 
in germanium so that failure of establishment 
of equilibrium in the electron populations of the 
valleys cannot be the source of the discrepancies. 
For example, the smallest intervalley scattering 
rate found in all of the samples studied by Wein- 
reich et al.° was 5x10® sec™ in a high-purity 
specimen. 

We thank Professor S. L. Quimby of Columbia 
University for permitting us to use certain 
components of his ultrasonic pulsing equipment. 





'We are indebted to Dr. G. E. Brock of the IBM 
Research Center for supervising the preparation of 
this unusually large and homogeneous specimen of 
heavily doped germanium. 

*H. J. McSkimin, J. Appl. Phys. 24, 988 (1953). 

*M. E. Fine, J. Appl. Phys. 24, 338 (1953). 

‘R. S. Leigh, Phil. Mag. 42, 139 (1951). 

5G. Weinreich, T. M. Sanders, Jr., and H. G. 
White, Phys. Rev. 114, 33 (1959). 
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EVIDENCE FOR THE DOUBLE F MODEL OF THE M CENTER 


Bruce J. Faraday, Herbert Rabin, and W. Dale Compton 


U.S. Naval Research Laboratory, Washington, D. C. 
(Received June 8, 1961) 


Recent results of van Doorn’ on the thermal 
equilibrium of F and M centers in KCl show that 
the concentration of M centers varies quadrati- 
cally with the concentration of F centers. This 
result supports the proposal of van Doorn and 
Haven? and Pick® that the model of the M center 
is a pair of associated F centers (F, model) rather 
than a vacancy pair plus a single F center as sug- 
gested earlier by Seitz* and modified by Knox.° 
This Letter presents new results on the relative 
concentration of F and M centers induced by ir- 
radiation, and the implication of these results on 
the M-center model. 

Thin plates of KBr and KCl, approximately 0.09- 
0.20 mm thick, have been irradiated at liquid he- 
lium and liquid nitrogen temperatures with x rays 
from a 50-kv de generator for periods up to sixty 
hours. Sensitive optical absorption measurements 
show that a weak M band appears with the genera- 
tion of F centers in KBr, KBr doped with Ca, KCl, 
and KCl doped with Sr. The data for KBr (Fig. 1) 
show that the concentration of M centers produced 
by these low-temperature irradiations is a quad- 
ratic function of the F-center concentration. A 
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FIG, 1. Absorption constant of the M band vs absorp- 
tion constant of the F band squared for x-rayed KBr at 
liquid helium and nitrogen temperatures. The optical 
measurements were taken at the irradiation tempera- 
ture. 


similar result was found in the KCl specimens. 

As is shown below, this is evidence that an M cen- 
ter is composed of two F centers. The fact that 
the divalent impurities and accompanying positive- 
ion vacancies did not alter the relative concentra- 
tions of F and M centers (see Fig. 1) is additional 
evidence that the positive-ion vacancy is not es- 
sential to M-center production. 

The x-ray production of M centers is predicta- 
ble from the probability of forming two F centers 
as neighbors. If it is assumed all F centers are 
produced with equal likelihood, the probability 
that two F centers will be formed as neighbors 
is proportional to (x;/N)’, where np is the con- 
centration of F centers and N is the concentration 
of negative-ion sites in the normal lattice. Thus 
the concentration of M centers is given by 

ny, =Kn"/N, (1) 
where K is a constant. Figure 1 shows this pre- 
dicted quadratic dependence. Using the Gaussian 
form of Smakula’s formula® for both F and M cen- 
ters, it follows that 


a WwW \/f.? K\{C_? 
(3! Wie M N Cu 


where C=y/(u?+2)?, uw being the index of refrac- 
tion, a =absorption constant in cm™, W=the ab- 
sorption half-width in ev, and f=the oscillator 
strength. The indices F and M refer to the F and 
M centers, respectively. Table I gives a tabula- 
tion of the half-widths and values of a,4/a;” ob- 
tained from a least-squares fit of the data of Fig. 1. 
Assuming that negligible diffusion of F centers oc- 
curs at liquid helium and nitrogen temperatures, 
M centers can only be formed during x raying 
when two F centers are formed as nearest neigh- 
bors. In this case K=12, corresponding to the 
number of nearest negative-ion sites surrounding 
an F center. Average values of f,, are taken from 
the literature.” The values of f,, so calculated 
(0.12 to 0.20, Table I) appear quite reasonable in 
view of the simplicity of the calculation.°® 

In light of the low-temperature results, two ex- 
amples of M-center formation at higher tempera- 
tures can be considered. First, the M-center con- 
centration in KBr and KCl crystals x-rayed at 
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Table I. Data for calculating f,), for KBr and KCl assuming K=12. 








2 a 
y/o p (cm) Wl" FE (ev/ev) LP tu 
KBr (He temp.) 1.87x10~ 0.052/0.158 3.2 0.15 
KBr (N, temp.) 1.71 x10 0.063/0.190 4.0 0.12 
KCl (He temp.) 1.951078 0.056/0.169 2.7 0.20 
KCl (N, temp.) 1.96 x10~% 0.066/0.194 2.9 0.18 





* average values of f, are taken from the literature for a Gaussian F band: fp ay(KBr) = 0.48 and (fp) ay(KCl) 


=0.53.8 


liquid helium or nitrogen temperatures increased 
by a factor of about fifty upon warming to room 
temperature.’° This suggests that F centers some 
distance apart aggregate to form M centers during 
warmup. Second, F- and M-center equilibrium 
after room-temperature x raying may be studied. 
The published data of Mador et al.” for NaCl and 
LiF are particularly appropriate for this purpose 
since optical absorption measurements were made 
transversely to the x-ray direction, insuring uni- 
form coloration along the optical path. A treat- 
ment of their data to test the quadratic relation 
between ay, and a, is shown in Fig. 2. It is again 
seen that a), is a linear function of a F° Equation 
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FIG, 2. Absorption constant of the M band vs absorp- 
tion constant of the F band squared for NaCl and LiF 
x-rayed and measured at room temperature, taken 
from Mador et al," 
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(2) can now be used to compute a value of K for 
NaCl and LiF at room temperature, assuming a 
value of the oscillator strengths.’* K may now be 
interpreted as follows. If two F centers are con- 
tained in a spherical volume of K negative-ion 
sites, they associate to form an M center. Ifa 
is the lattice constant, i.e., the distance between 
like ions along a (100) direction, then a “sphere 
of interaction” for F centers of radius 7 is given 
approximately by the relation (417 /3)7>~Ka*/4. It 
follows that r/a™0.39(K)”°, which is the maximum 
number of lattice constants over which F centers 
join to form M centers. Using values of K calcu- 
lated from Eq. (2) from the least-squares slopes 
of Fig. 2, r/a can thus be determined. Values of 
K and r/a are as follows: 





K r/a 
NaCl 7000 7.4 
LiF 4400 6.4 


This result implies that the interaction of two F 
centers is sufficiently great at room temperature 
that two F centers separated by six or seven lat- 
tice constants join to form M centers. The nature 
of the mechanism which promotes such an associ- 
ation is not understood. Alternatively this might 
be taken as evidence for a nonrandom distribution 
of F centers. 

A value of K for KCl can also be determined for 
van Doorn’s high-temperature data’ (~ 1000°K). 
This value is roughly a factor of ten lower than 
the room-temperature values, suggesting an equi- 
librium between the formation and the breakup of 
M centers at the high temperature. 

A substantial amount of evidence now appears 
to support the F, model of the M center. It has 
been demonstrated that F and M centers are quad- 
ratically related in cases of both thermal quench- 
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ing and x raying over a wide temperature and con- 
centration range in a variety of alkali halides. 
Furthermore, this model is in accord with the 
requirement of inversion symmetry,** the absence 
of dielectric loss'* and magnetic susceptibility,”® 
and the inconclusiveness of electron spin reso- 
nance’® of the M center. 

In light of the above discussion the suggestion 
of Pick® that the R and N centers are higher ag- 
gregates of F centers has been examined. Treat- 
ing unpublished data’’ for NaCl that had been uni- 
formly irradiated with gamma rays at room tem- 
perature, it was found that the R,-center absorp- 
tion varied approximately as the 1.5 power of the 
M-center absorption. This supports the conten- 
tion that the R center is an aggregate of an M and 
F center, or simply three F centers as suggested 
by Pick. A more extensive study is indicated. 

The authors are indebted to Dr. C. C. Klick for 
many helpful discussions. 
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SIZE-DEPENDENT RELAXATION TIMES AT LIQUID HELIUM TEMPERATURES* 


F. R. Nash 


Columbia Radiation Laboratory, Columbia University, New York, New York 
(Received June 26, 1961) 


Recovery times have been observed in one of 
the copper Tutton salts which depend upon the 
size of the crystal(s) employed.’ The salt in 
question is fully concentrated Cu(NH,),(SO,),- 6H,O. 
The paramagnetic ion is Cu*® (S = ) and the two 
ions per unit cell give rise to two lines. How- 
ever, at high concentrations the two lines coa- 
lesce to form one broad line whose full width at 
half maximum varies from 170 to 450 oersteds 
depending upon the angle between the crystal 
axis and the magnetic field. 

The relaxation times were measured with an 
X-band pulse-decay superheterodyne spectrom- 
eter, the output of which was the dynamic re- 
covery curve of the spin system.* The relaxa- 
tion times were determined from the exponential 
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tails of the decays, in the region where the spin 
temperature = lattice temperature. 

The crystals reported on were prepared in two 
different ways. In the first instance a fine pow- 
der of crystallites was allowed to grow slowly 
over a period of 24 hours. The second method 
involved the preparation of a supersaturated so- 
lution by heating. The solution was stirred con- 
tinually as it cooled to room temperature with 
the result that a fine powder of tiny crystals was 
deposited in the bottom of the beaker. In both 
cases the crystals were then dried by exposure 
to air and carefully sorted according to size 
using a series of sieves. The hole sizes in the 
ten sieves ranged from 37 to 1000 microns. 

Only naturally grown single crystals were used 
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in these experiments although some single crys- 
tals were crushed and tested to see that the size 
effect was qualitatively unaltered. 

The cavity was entirely filled with Styrofoam 
except for a pocket = 10 mm* which was located 
in the region of the cavity where H,,=0. This 
placement was to allow the use of an appreciable 
volume of sample and still keep the cavity loading 
<10%.° 

The results are presented in Figs. 1 and 2. It 
seems reasonable to suggest that crystals which 
are grown quickly (=1 minute) are sufficiently 
strained so that some internal energy conver- 
sion mechanism alleviates the apparent bottle- 
neck. One might suggest phonon-phonon inter - 
actions but the observed temperature dependence 
is not as strong as this mechanism would imply. 

Within the experimental error of =~+10% one 
can see that T; /,<L, where L is some charac- 
teristic length of any one of the many crystal- 
lites in a given portion of sample. 

From measurements made on single crystals 
whose smallest dimensions were = 1-2 mm, it 
was noted that the transition from the region 
where 74 /e «L to the region where 7 /e is in- 
dependent of LZ occurs for L=1 mm. Further- 
more, the relaxation times of the single crystals 
were =20% shorter than those times in the sam- 
ple where L=855 microns. Presumably this is 
related to the fact that large single crystals are 
somewhat more strained than smaller ones. The 
same kind of bottleneck has been seen in the 
deuterated crystals of the same species and 
concentration and also in 30% concentration 
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FIG. 1. Relaxation times in quickly grown crystals 
of 100% Cu(NH,)2(SO,).°6H,O vs characteristic length 
of crystal. 
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hydrated crystals. 

In order to test that the liquid helium in con- 
tact with the sample was really cooling the in- 
dividual particles, a very light smear of the L 
<37 microns sample was made on the surface 
of the Styrofoam cavity insert in the region of 
(Hrf)max: This was to assure the helium liquid 
clear access to all the grains in the sample. 
Exactly the same results were obtained as when 
the sample consisted of 10 mm‘ of tightly packed 
powder, in the weak-field region of the cavity. 

Spatial spin diffusion as the explanation for 
this size effect can be ruled out on the basis 
of an order-of-magnitude calculation of the over- 
all diffusion time as well as dependence upon 
length and lattice temperature.* The theories 
of phonon imprisonment in paramagnetic crys- 
tals® would indicate that lattice-bath relaxation 
times are being observed.® Indeed, there is 
fairly good agreement between the theory® and 
the observed results. To see that phonon trap- 
ping is reasonable one can make an order -of- 
magnitude estimate of the phonon mean free 
path (m.f.p.) in the central region of the spin 
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FIG, 2. Relaxation times in slowly grown crystals 
of 100% Cu(NH,)2(SO,).°6H,O vs characteristic length 
of crystal. 
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resonance: 


m.f.p. = (nm) yePv)avT, /AN 


= (kT /hv)?240V*AvT,, /v°N, 


where (n) ay =number of phonons per mode =kT/ 
hv, v is the mean acoustic velocity, p(v) is the 
density of lattice modes =127v*/v*, Av is the 
full width of the spin resonance, T, is the spin- 
lattice relaxation time, and AN=NkT/2hv. If 
T=1.34°K, v=9.6x10° cps, T,=10~ sec, v =2.5 
x10° cm/sec, and N =3 x10" spins/cm’, then 
m.f.p.=30 microns or 3x10 cm. 

There are, however, two aspects of this effect 


which are not clear. 
(1) Tests were made on a fully concentrated 


single crystal in which the sample was cooled 
by helium gas (pressure=1 mm of Hg) at 4.21°K, 
2.09°K, and 1.34°K. The run was repeated but 
this time the crystal was cooled directly by liq- 
uid helium. Within experimental error the re- 
laxation times in both cases were the same. The 
test was performed on another single crystal of 
a different shape and again it did not matter 
whether the crystal was cooled by helium liquid 
or gas. 

(2) Phonon trapping theoretically® depends upon 
a direct spin-lattice interaction. In addition, ex- 
periments which have been performed by the 
author at lower concentrations indicate that 
the spin-lattice relaxation time is concentration 
dependent. However, the theoretical calcula- 
tions’ which have been carried out for the con- 





centration-dependent, first order spin-lattice 
interactions are inadequate to explain the short- 
ness of the spin-lattice relaxation times obtained 
by extrapolating the present data to L=0. 

The author would like to acknowledge his in- 
debtedness to Professor C. H. Townes and to 
Dr. J. A. Giordmaine for their innumerable sug- 
gestions, their time, and their encouragement. 

Thanks are also due to Professor H. Boorse 
and H. Louie for liquid helium and to B. Biavati 
and S. Zemon for discussions concerning low- 
temperature techniques. 
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Office of Scientific Research. 
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OSCILLATORY PHOTOCONDUCTIVITY IN InSbt 


W. Engeler,* H. Levinstein, and C. Stannard, Jr. 
Syracuse University, Syracuse, New York 
(Received June 21, 1961) 


The addition of Cu to InSb has been shown by 
Blunt! to produce two deep acceptor levels. He 
also observed a series of oscillations in the ex- 
trinsic photoconductive response associated with 
the lower of these two levels. To investigate the 
mechanism responsible for this behavior the prop- 
erties of impurities similar to Cu were examined. 
Ag and Au, which had also been found to produce 
two deep acceptor levels,? were examined in ad- 
dition to Cu and Ni. 

Oscillatory behavior of the photoconductive re- 
sponse at liquid He temperature was observed 
for all of these impurities as shown in Fig. 1. In 
all cases the oscillations observed have minima 
separated by constant energy intervals. The loca- 
tions of these minima depend on the impurity and 
are thus distinct from lattice absorptions at 27.5 
and 33.5 uw which are independent of impurity. 
Figure 2 shows a plot of the number of the mini- 
mum vs the energy at which it occurs, for all of 
the impurities investigated. A straight line is ob- 
tained for each impurity. The slopes of these 
lines, corresponding to the energy difference be- 
tween successive minima as calculated by a least- 
squares fit of the data, are identical within ex- 
perimental error. In the case of gold it was pos- 
sible to obtain the value of the long-wavelength 
threshold of the photoconductive response. The 
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FIG, 1. Extrinsic photoresponse of InSb. 
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intercept of the straight line with the energy axis 
yields an energy value in agreement with the 
threshold wavelength. This value is indicated by 
the arrow in Fig. 1. Since the threshold of the 
Cu and Ag response curves was not within the 
range of the optical equipment available, the as- 
signment of the number of the minima was based 
on the estimate of the position of threshold ener- 
gy by Hall measurements. The fourth curve shown 
in Fig. 1 was obtained from a sample which was 
originally prepared to contain Ni. The photocon- 
ductive response cutoff of this sample shows a 
wide departure from the energy of the Ni level as 
determined by the Hall data. It is highly probable 
therefore that the response is due to an unknown 
impurity which was introduced into the sample 
during its processing. Since the maximum solu- 
bility of Ni in InSb is of the order of 10** atoms/ce, 
only a small additional impurity would be neces- 
sary to cause such an effect. The values of slopes 
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purities in InSb. 
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Table I. Summary of results. 
Activation energy (ev) 
Curve Impurity Slope (ev)# Intercept (ev) (from Hall measurements) 
1 Cu 0.0244 0.026 0.023 
2 Ag 0.0244 0.028 0.023 
3 Au 0.0244 0.043 0.032 
4 0.0244 0.035 ose 





*prror does not exceed 0.0001. 


and intercepts together with Hall data for the vari- 
ous impurities are summarized in Table I. 

These data strongly suggest that the observed 
spectral response is composed of a series of 
superimposed spectral curves. The long-wave- 
length threshold corresponds to the excitation of 
a hole from the impurity state to the valence band. 
The next threshold corresponds to the same proc- 
ess accompanied by the emission of a longitudinal 
optical phonon, succeeding thresholds at shorter 
wavelengths to the emission of 2, 3, 4, etc., 
phonons, respectively. 

The slopes of the curves in Fig. 2 are then the 
energy of the longitudinal optical phonon. The 
value obtained in this manner is in good agree- 
ment with the value of 0.024 ev obtained by Hall, 
Racette, and Ehrenreich® from tunnelling experi- 
ments and the value of 0.0244 ev recently reported 
by Fray, Johnson, and Jones.* 

A process of this type can be explained as fol- 
lows. Starting from the threshold at k=0 and 
going to higher energy, the transition probability 
for direct transitions to the impurity level first 
increases with the increase in the density of states 
of the valence band. This is followed by a decrease 
of this transition probability for larger k values 
due to the decrease in the density of states of the 
impurity level. This decrease is slowed by the 
possibility of indirect transitions. When a photon 
energy is reached where a longitudinal optical 


phonon can be emitted, the transitions can again 
take place at k=0 and the process is repeated 

with the emission of this phonon. The optical ab- 
sorption constant is directly related to this transi- 
tion probability. The magnitude of the spectral 
response in turn depends on the relative magni- 
tude of this absorption and other competitive ab- 
sorption processes. Because of the lack of knowl- 
edge of the competitive processes, the magnitude 
of the observed oscillations cannot be directly re- 
lated to the transition probabilities. This behavior 
should be observable in other polar compounds and 
should be a valuable tool in the study of these ma- 
terials. 

The authors acknowledge the assistance of L. 
Apker, H. Ehrenreich, and W. W. Tyler of the 
General Electric Research Laboratory in the in- 
terpretation of these results. 





+tSponsored by Air Research and Development Com- 
mand, U. S. Air Force. 

*Now at General Electric Research Laboratory, 
Schenectady, New York. 
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*w. Engeler and H. Levinstein, Bull. Am. Phys. 
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5R. N. Hall, J. H. Racette, and H. Ehrenreich, 
Phys. Rev. Letters 4, 456 (1960). 

4s, J. Fray, F. A. Johnson, and R. H. Jones, Proc. 
Phys. Soc. (London) 76, 939 (1960). 


63 








VOLUME 7, NUMBER 2 


PHYSICAL REVIEW LETTERS 


Jury 15, 1961 





AUTO-IONIZATION OF THE (1s2s2p) *P,, LEVEL IN Li AND He 


Jerry L. Pietenpol* 


Columbia University, New York, New York 
(Received June 20, 1961) 


The negative helium ion has been observed ex- 
perimentally,*-* and Sweetman‘ has obtained a 
lower bound of 10~ second for its lifetime. 
Arguments have been advanced®”® indicating 
that the observed metastable state is a (1s2s2p) 
“Pes configuration. The purpose of this note is 
to report a theoretical estimate of the lifetime 
of the state for auto-ionization. The calculation 
has been made both for negative helium and for 
lithium; the result is of interest in the latter 
case because of the possibility of measuring 
the fine and hyperfine structure of the level. 

The state is radiatively metastable because 
of the spin selection rule. Holgien and Midtdal® 
have shown by an elaborate variational calcula- 
tion that this level in He” lies below the 1s2s 
level of the helium atom, and have pointed out 
that the state is therefore not subject to auto- 
ionization in the usual way. The only available 
final state consistent with angular-momentum 
and parity conservation is then (1s*kf) *F,.. The 
Coulomb, spin-orbit, and spin-other-orbit oper- 
ators have vanishing matrix elements for the 
transition, but the process can proceed by the 
magnetic spin-spin interaction between two elec- 
trons. 

If we write the spin-spin operator in the form 


H' = -Uo'(6,* V,) (G2 V,)(1/%42), 


the matrix element can be evaluated by a double 
partial integration. (It is necessary to exclude 
an infinitesimal sphere about the point of singu- 
larity, so that one must look for surface terms 
here. These terms vanish, however, for any 
matrix element between states of different spin.) 
We used antisymmetrized product wave functions, 
with screened hydrogenic functions for the bound 
electrons and a plane wave for the continuum. 
The effective Z’s and the energy of the initial 
state were obtained from a variational calcula- 
tion, including exchange between 2s and 2p. The 


parameters are as follows: 

For Li: initial-state Z,,.=3.0, Z9,=2.0, Zp 
=1.7, final-state Z,, =2.69, kinetic energy of 
continuum electron E =3.95 rydbergs. 

For He’: initial-state Z},=2.0, Z9,=1.0, 
Z2p = 0.67, final-state Z;,.=1.69, E=1.45 ryd- 
bergs. 

The integrations were carried out analytically. 
The result for the lifetime of the state is 1.6 
x107° sec for Li, and 1.7x10~° sec for He™. 
Variation of the Z parameters within reason- 
able limits produces a change of no more than 
a factor of two in the lifetime. In the very sim- 
ple approximate wave functions employed, the 
use of a plane wave for the continuum state is 
justified for helium (since the centrifugal bar- 
rier effectively prevents the electron from pen- 
etrating the 1s” core), but is questionable for 
lithium. The initial-state wave function did not 
include correlation between the 2s and 2p elec- 
trons; this effect would be expected to lengthen 
the calculated lifetimes. 

The author is indebted to Professor H. M. Fo- 
ley and Professor R. Novick for valuable discus- 
sions, and to Professor N. M. Kroll for suggest- 
ing the spin-spin interaction as the mechanism 
of the process. 
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792 (1956) [translation: Soviet Phys.-JETP 3, 764 
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STRUCTURE IN THE KINETIC ENERGY SPECTRUM OF FRAGMENTS FROM 
THERMAL-NEUTRON-INDUCED FISSION OF U3*t 


W. M. Gibson 
Bell Telephone Laboratories, Murray Hill, New Jersey 


and 


T. D. Thomas and G. L. Miller 
Brookhaven National Laboratory, Upton, New York 
(Received June 19, 1961) 


We have observed fine structure in the kinetic 
energy spectra of the heavy fragments from ther- 
mal-neutron-induced fission of U**. This fine 
structure in the energy spectra is reflected as 
fine structure in the calculated fission-mass- 
yield distribution. 

A target of 5 ug/cm? of U**® evaporated onto a 
15-.g/cm? backing of VYNS film was placed be- 
tween two 1.5-cm? p-n junction semiconductor 
particle detectors 5 cm apart. The target, coun- 
ters, and integrating preamplifiers were enclosed 
in an evacuated aluminum can, which was lowered 
into the thermal column of the Brookhaven graphite 
research reactor. The pulses from the preampli- 
fiers were further amplified by linear amplifiers 
whose outputs were fed into a coincidence circuit 
and into the two axes of a 64 x64 two-dimensional 
pulse-height analyzer.’ The pulse-height analyzer 
was triggered by the coincidence circuit and, there- 
fore, recorded only those events in which a fission 
fragment entered each counter within the resolv- 
ing time of the coincidence circuit (~0.1 psec). 

Two separate experiments were performed. In 
the first experiment, counters made by diffusion 
of phosphorus into 100-ohm-cm p-type silicon 
were used at a reverse bias of 100 volts. The 
coincidence counting rate was 2400 fissions/min- 
ute, and a total of 1.66 10° fissions was recorded. 
In the second experiment, the counters were of 
300-ohm-cm silicon biased at 100 volts. Because 
there was evidence for a small low-energy tail on 
the kinetic energy spectra in the first experiment, 
two Plexiglas collimators were used in the second 
to keep the fission fragments from hitting the edges 
of the counters. The coincidence counting rate was 
500 fissions/minute with a total of 6.1 x10° fissions. 
Summation of the mass spectra over the total kinet- 
ic energy range of 70 Mev to 250 Mev gave a peak- 
to-valley ratio of 370/1. 

Figure 1 shows the very pronounced structure 
in the energy spectra of heavy fragments in co- 
incidence with light fragments of greater than 100 
Mev. The energy channels are about 1 Mev wide 
in both directions. Stein’s time-of-flight data have 


been used to calibrate the energy scale.? Recent 
time-of-flight measurements on thermal-neutron 
fission of U*** by Milton® show essentially the same 
structure. We observed no structure in the spec- 
tra of the light fragments nor in the spectra of 
heavy fragments in coincidence with light frag- 
ments of less than about 100 Mev. 

The fission yield and average total kinetic energy 
as a function of mass number calculated from the 


4 
10 T “T T T sa T T 
5 





105 


{ 106 £ 
| [i ( \ 
| ii ij } . ; 
i i 7 } AX 

slid tJ Ale 


HEAVY FRAGMENT ENERGY IN Mev 


we T 
— 
ee 
~~ 


~—— 
i 


COUNTS PER CHANNEL 
ro) 
w 

howl i 0 | 


i 








ny 





FIG. 1. Spectra of heavy fission fragments in coin- 
cidence with given light fragment energies, E;. The 
open points and broken curves correspond to a heavy 
fragment into counter A and a light fragment into 
counter B; the closed points and unbroken curves to 
the reverse. 


65 








VOLUME 7, NUMBER 2 





PHYSICAL REVIEW LETTERS 


Jury 15, 1961 

















3 
180 & 
Wd 
on gan & 
r ° ° °% 170 | 
~ e e *%, < 
ra 3 *m, — 11600 
“ws” A se “ey, a 
we “J 150 9 
e ° x 
» A ~ 
6} 4140 
= < 
a 
Oo St : 
a 
Ww 
ag 4 
z | 
93 
” 
> ah : 
z 
° 
W 1b 4 
vn 1 
we 
16 80 90 100 110 120 130 140 150 160 
MASS NUMBER 
FIG, 2. Fission yield and average total kinetic en- 


ergy as a function of mass number. 


experimental data are shown in Fig. 2. Except in 
the regions of mass 100 and mass 134, the mass- 
yield data are in good agreement with the data ob- 
tained by Milton, who used the time-of-flight tech- 
nique.* The kinetic energy curve falls somewhat 
lower than does Milton’s, the difference being due 
to a difference between the time-of-flight meas- 
urements of Milton and those of Stein.” 

The two shoulders on the mass-yield curve oc- 
cur at light fragment masses of 92 and 101 and 
the peak occurs at light mass 96; the maxima in 
the curves of Fig. 1 correspond to light fragment 
masses of 91, 96, and 101. Thus, the structure 
in the mass-yield curve appears to be related to 
the structure in the kinetic energy spectrum. 
Mass-yield curves for specific values of the total 





kinetic energy differ from the summed curve of 
Fig. 2 in that they are narrower and in general 

do not show the fine structure apparent in Fig. 2. 
The structure in the total mass distribution arises 
from the fact that the mass-yield curves for dif- 
ferent total kinetic energy tend to fall into three 
groups. For total kinetic energies less than 162 
Mev, the most probable light fragment mass num- 
ber is about 92; for energies between 164 and 171 
Mev, about 96; and for energies greater than 173 
Mev, about 101. 

The highest energy peaks in the curves shown 
in Fig. 1 occur at energies corresponding to a 
heavy mass of about 135, which may be associated 
with the closed shells of 82 neutrons and 50 protons. 
However, there is no obvious correlation between 
closed shells and the other two peaks which cor- 
respond to heavy fragment masses of about 140 
and 145 (or light masses of 96 and 91). 

We have conducted similar measurements on 
the spontaneous fission of Cf**. An experiment 
in which 6 x10° fissions were recorded showed 
that any fine structure that may exist in the kinetic 
energy spectrum of Cf*™ is much less pronounced 
than in the thermal neutron fission of U***. 

We are pleased to acknowledge the assistance of 
C. A. Boulin, who constructed the preamplifiers, 
and T. C. Madden, who made the counters and 
helped design the apparatus. 





TResearch was performed in part under the aus- 
pices of the U. S. Atomic Energy Commission. 

'R, L. Chase, Brookhaven National Laboratory Re- 
port BNL-3838, 1958 (unpublished). 

*w. E. Stein, Phys. Rev. 108, 94 (1957). 

33. C. D. Milton and J. S. Fraser, following Letter 
[Phys. Rev. Letters 7, 67 (1961). 

43. C. D. Milton, quoted by G. C. Hanna, Proceed- 
ings of the International Conference on Nuclear Struc- 
ture, Kingston, Canada, edited by D. A. Bromley and 
E. W. Vogt (University of Toronto Press, Toronto, 
1960), p. 866. 
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PROMPT FISSION YIELDS AND TOTAL KINETIC ENERGY BEHAVIOR 
FROM TIME-OF-FLIGHT MEASUREMENTS 


J. C. D. Milton* 
Lawrence Radiation Laboratory, Berkeley, California 


and 


J. S. Fraser 


Chalk River Laboratories, Atomic Energy of Canada Limited, Chalk River, Ontario, Canada 
(Received June 19, 1961) 


An improved version of the time-of-flight appa- 
ratus described previously’*? has been used to 
study the thermal neutron fission of U***, U%, 
and Pu**® using a neutron beam of intensity about 
10° neutrons sec™ cm™ from the NRU reactor. 
The sources, ranging from 20 to 40 ug cm™, 
were electrosprayed onto two cross-laminated 
VYNS films of total surface density 7 ug cm~ 
which had been previously coated with 18 ug cm~? 
of gold to make them electrically conducting. The 
energy loss caused by source and backing thick- 
ness was measured for the actual sources used 
in the experiment and found to be 1.1+ 0.3 Mev 
in good agreement with the value calculated from 
the known thicknesses. 

The fragments were detected at the end of their 
nearly equal 178-cm flight paths by 4-inch diam- 


eter secondary electron detectors which gave the 
apparatus a resolving time (FWHM) of 2.5 nano- 
seconds. The mass resolution including the con- 
tribution of timing jitter, recoil from neutrons, 
and source thickness was calculated to be 2 mass 
units. 

In all, well over 10° events were processed on 
U**> +n (including some with the unequal flight 
paths of 270-180 cm) but the most carefully con- 
trolled series consisted of 1.48x10°, 1.55 x10°, 
and 1.03 x10° events for U**, U*°5, and Pu”, 
respectively, and it is on this set that the data 
presented here are based. The experimental 
results on Cf”? are not new but are those of 
reference 1. 

An experiment of this type conveys an enormous 
amount of information. The purpose of this Let- 
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FIG, 2. Mass-energy contours for U**, U5, and 
Pu”**, The numbers attached to the contours are events 
(10° fissions)~! (mass unit)~! (2.5 Mev)~!. Some idea 
of the statistical accuracy may be gained from the fact 
that the number of events in each cell of the matrix 
from which the contour diagram was made is approx- 
imately 24 times the contour label. 
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ter is to note 3 points: (1) There is a dramatic 
drop in the total kinetic energy near symmetry; 
(2) there is a clearly marked fine structure in 
the mass yield which is most prominent at high 
total kinetic energies; and (3) a possible explana- 
tion for the phenomenon (1). 

The behavior of the kinetic energy is illustrated 
in Fig. 1. Except in the case of Cf*™, there is a 
drop in kinetic energy of 30-40 Mev in approach- 
ing symmetric mass division. At first it seemed 
impossible to account for this dip in terms of the 
balance of energy between the mass energy re- 
lease E R? and the excitation energy made up of 
components Ep (neutron binding), E,, (neutron 
kinetic), and E, (gamma rays), but this turned 
out to be incorrect. The solid points in Fig. 1 
are calculated from the formula Ex =max(Ep) 
-v((Ep) +E pp) -E,. The quantities E,, and Ey 
were assumed to be independent of mass and giv- 
en the approximate values E,, p® 1.2 Mev, E, = 7.5 
Mev. The averages were taken over a Gaussian 
charge distribution with an rms deviation of 0.7 
charge unit. Using the mass formula of Cameron,’ 
(Ep) was found at the point where (Ep) had its 
maximum value. These quantities have been cal- 
culated for 25 fissioning species between Ra” 
and Fm™ using an IBM-709 computer. However, 
the principal difference between the E K calcula- 
ted here and those done previously lies in the 
selection of v. It was assumed that v depends 
only on the mass of the fragment and not on the 
fissioning nuclide. This may be seen to be a good 
assumption by looking at the values of v calculated 
by Terrell.* Accordingly, a smooth curve was 
drawn through Terrell’s family of curves. The 
values of ) given by this smooth curve are 2.50, 
2.63, 2.86, and 3.81 for U**, U*5, Pu?%, and 
Cf*™, respectively. This procedure makes v 
available for all masses between 82 and 159 with 
the exception of the region 120 to 130. In this re- 
gion the two ends of Terrell’s curve were joined 
by a straight line. At first sight this may seem 
to be a very shaky business but, in fact, almost 
any interpolation between Terrell’s two curves 
will do since in U*** (for example) it affects the 
shape of the kinetic energy between mass 130 and 
118 but not the value at these two masses. There 
would still be a drop of nearly 40 Mev, even if v 
were taken as identically zero for the region 121- 
130. 

The total excitation energies at symmetry pre- 
dicted by this picture are large, about 48 Mev for 
U*** and about 6 neutrons would be expected. The 
existence of two types of fission—symmetric fis- 
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sion with highly excited fragments of low kinetic 
energy and asymmetric fission with moderately 
excited fragments of high kinetic energy— might 
be associated with the existence of two barriers 
and two saddle-point configurations in the liquid 
drop calculations of Cohen and Swiatecki.® The 
height of the second barrier relative to the first 
falls rapidly as x increases above a critical value 
in the vicinity of 0.7 and its relatively low energy 
and apparent stability against asymmetry® suggest 
identifying this barrier with symmetric fission. 
All three of the fissioning species showed fine 
structure in the mass yield between masses 134 
and 146. This fine structure shows up in the total 
yields® but it is perhaps more instructive to look 
at the contour diagram, Fig. 2. Here it is clearly 
seen that the structure appears predominantly at 
high kinetic energies.’ It seems to be more a re- 
sult of inhibition than of preference and must be 
closely related to the fact that the events in this 
region are “running out of energy” so that the 
nuclei formed are nearly in their ground states. 
The contours along the northeast face of the moun- 
tain are nearly parallel to the line of maximum 
energy release. The structure is less pronounced 


in Pu**® and scarcely discernible in Cf***. In these 
cases the mountain is further away from the lim- 
iting energy. 

A full account of this work will shortly be sub- 
mitted to the Canadian Journal of Physics. I am 
indebted to Dr. W. J. Swiatecki for pointing out 
the connection with the liquid drop calculations. 








*On leave of absence from Atomic Energy of Canada 
Limited, Chalk River, Ontario, Canada. 

‘J. C. D. Milton and J. S. Fraser, Phys. Rev. 111, 
877 (1958). 

23. S. Fraser and J. C. D. Milton, Nuclear Instr. 2, 
275 (1958). 

3a. G. W. Cameron, Can J. Phys. 35, 1021 (1957). 

‘James Terrell, Bull. Am. Phys. Soc. 6, 16 (1961); 
and private communication. 

5§. Cohen and W. J. Swiatecki, Aarhus University 
Report, January, 1961 (unpublished). 

83. C. D. Milton, Proceedings of the International 
Conference on Nuclear Structure, Kingston, Canada, 
edited by D. A. Bromley and E. W. Vogt (University 
of Toronto Press, Toronto, 1960), p. 866. 

"This is also shown in the work of W. M. Gibson, 

T. D. Thomas, and G. L. Miller, preceding Letter 
[Phys. Rev. Letters 7,65 (1961)]. 











DEUTERONS FROM HIGH-ENERGY PROTON BOMBARDMENT OF MATTER 


S. T. Butler and C. A. Pearson 
The Daily Telegraph Theoretical Department, School of Physics, University of Sydney, Sydney, Australia 
(Received June 16, 1961) 


The mass analysis of the secondary particles 
from targets exposed to the 25-Bev proton beam 
of the CERN proton synchrotron has revealed an 
appreciable fraction of high-energy deuterons.’ 

In addition large numbers of deuterons have been 
reported among the “grey” tracks in nuclear emul- 
sions.?* 

Recently Hagedorn® has suggested that these 
deuterons are produced in elementary nucleon- 
nucleon collisions. However, the results obtained 
from a statistical model are strongly energy de- 
pendent and especially at lower energies it is hard 
to reconcile them with the experimental results. 
For deuterons with momenta below 2 Bev/c they 
disagree with experiment by a factor greater than 
10°, 

In interpreting the emulsion results the Bristol 
group suggested that the deuterons are formed 
when a knock-on nucleon pairs with a nucleon in 
the tail of the Fermi distribution of the excited 
Nucleus. Because of their high energy, however, 
itseems more likely that they result from the 
pairing of two cascade nucleons. The angular 


distribution of the knock-on “shower” nucleons 
has a strong forward peak and in many cases 
there will be pairs of particles which have small 
relative momenta. It is from such pairs that we 
believe deuterons are formed. 

To calculate the ratio of deuterons to protons 
expected from such a process, we have considered 
the following model. An incident beam of free 
protons and neutrons is scattered by an effective 
nuclear potential, which represents the interaction 
of the high-energy cascade nucleons with the re- 
mainder of the nucleus. The protons and neutrons 
can interact with each other during the scattering. 
The incident particles are not parallel or mono- 
energetic, but have a certain momentum distri- 
bution, say P(k)dk, which is the same for neutrons 
and protons and fits the experimental distribution 
for cascade protons. The flux of particles in the 
incident beam is chosen to yield n, nucleons of 
each type contained in the nuclear volume at any 
one time, where n, is the experimental average 
number of high-energy protons from a struck 
nucleus. 
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FIG. 1. Diagrams (a), (b), and (c) of Fig. 1 illus- 
trate the simplest methods of deuteron formation. ky, k, 
are the momenta of the proton and neutron in the initial 
state, q the recoil of the nucleus, and K the deuteron 
momentum in the final state. In case (a) the proton is 
scattered into an intermediate state by an interaction 
with the nucleus. The scattered proton and an unscat- 
tered neutron then interact with each other to form a 
deuteron. In case (b) a scattered neutron pairs with an 
unscattered proton. In case (c) the neutron and proton 
interact first with each other to form an intermediate 
deuteron state. This deuteron is then scattered by the 
nucleus into the final state. 


To form a deuteron, a proton and a neutron must 
interact with each other and with the nuclear field, 
so that in terms of Feynman-like diagrams the 
simplest schemes of deuteron production can be 
illustrated as in Fig. 1. Here k,,k, represent 
the momenta of the proton and neutron, K the 
deuteron momentum, and q the recoil of the nuc- 
leus. In diagrams (a) and (b) one of the nucleons 
has been scattered before the two interact to form 
a deuteron; in (c) they interact first with each 
other and the virtual deuteron formed is scattered 
into the final state. 

The contribution from these diagrams has been 
calculated using a relativistic second order per- 
turbation theory developed for this problem which 
will be described in a later paper. In the rela- 
tivistic region the free-particle states are repre- 
sented by the positive-energy solutions of the 


Klein-Gordon equation and the nuclear field by 
the effective potential calculated by Duerr® for 
such particles. In the nonrelativistic region the 
nucleus is represented by the well-known optical 
potential. 

Since the formation of deuterons arises only 
from pairs of particles with small relative mo- 
menta, the interaction between the neutron and 
proton can be represented in their c.m. frame 
by the ordinary deuteron potential. 

To simplify the calculation we have chosen K 
to be in the forward direction, and for comparison 
with the experimental results we have calculated 
the ratio of the number of deuterons expected ina 
solid angle dQ about K with the number of free 
protons of momentum K which have passed through 
the nuclear potential and which emerge in the same 
solid angle dQ. 

In the calculation the following parametric form 
for the momentum distribution P(k) was chosen: 


P(k) ak (cose)”, (1) 


where a and § are parameters and @ is the angle 
of the vector k with respect to the direction of the 


incident nucleon which initiates the nucleon shower. 


This momentum distribution is normalized such 
that the average number of nucleons of each type 
within the nuclear volume at any one time, and of 
momenta exceeding a minimum value k,, is the 
number n,. In the CERN experiment, for example, 
the high-energy “cascade” protons have energies 
mostly exceeding 500 Mev, and their momentum 
distribution above this energy can be fitted by a 
distribution of type (1); in this case therefore we 
choose k, to correspond to an energy of 500 Mev. 

In general, from an accumulation of experimen- 
tal emulsion observations*»’~*° the parameter val- 
ues listed in Table I were chosen. The values of 
n, in this table refer to the silver and bromine 
nuclei of emulsions and would be expected to vary 
somewhat with the size of the nucleus. 


Table I. Values of the parameters a, £, and my for different minimum energies in the nucleon cascade. These 


values have been chosen from an accumulation of experimental emulsion observations, 


a 





Range of deuteron energies calculated 





Minimum energy (#*k,?/2m) in the nucleon cascade a B no 
30 Mev 3.5 1 6 60-250 Mev 
125 Mev 3.5 2 4 250-500 Mev 
250 Mev 4 3 2 >500 Mev 





*See references 4, 7-10. 
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FIG, 2. A comparison of the observed and calculated 


deuteron to proton ratio for protons and deuterons of the 
same momentum. Curve 2 represents the experimental 
results from emulsions.*-* Curves 3 and 4 are the deu- 
teron to proton ratios found in the CERN experiment. 
Curves 5 and 6 are the results calculated by Hagedorn 
for the two different deuteron factors Qg=0.2, 0.1, re- 
spectively. Curve 1 represents the deuteron to proton 
ratio calculated from the model described above. 


The results of this calculation are plotted in Fig. 
2 and compared with the experimental results of 
emulsion groups and the CERN group. The im- 
portant conclusion is that there are ample numbers 
of deuterons formed by this process to account for 
the experimental observations. Indeed, in the 
high-energy region the numbers yielded by the 
present calculation are considerably too large; 
this is almost certainly due to the large repulsive 
Duerr potential® which we have simply taken at 
face value. The deuteron ratio is proportional to 
the square of the modulus of the effective nuclear 
potential, and is therefore quite sensitive to this 
potential. It would in fact be possible that the 
present deuteron formation mechanism can be 
used to provide a direct and sensitive measure - 


ment of the strength of the effective nuclear po- 
tential for relativistic particles. 

The results of Hagedorn’ are also shown in Fig. 
2 for comparison purposes. It is difficult to rec- 
oncile the strong energy dependence of this model 
with the rather slow energy variation of the ex- 
perimental results. Moreover, it appears to the 
present authors extremely probable that the sta- 
tistical model employed by Hagedorn would 
strongly overestimate the deuterons produced 
in elementary nucleon-nucleon collisions at all 
energies, thereby impairing the “local” high- 
energy agreement of Hagedorn results with ex- 
periment. 

In the CERN experiment moreover the percent- 
age of deuterons detected increases with the size 
of the target nucleus roughly as the nuclear radius. 
If the deuterons are produced in nucleon-nucleon 
collisions, this is difficult to explain. On the pres- 
ent model involving the recombinations of cas- 
cade nucleons themselves such an increase of the 
deuteron ratio with the size of the target nucleus 
is to be expected. 

Full details of the present calculation will be 
described in a later paper. 

The authors are grateful to Professor Messel 
for providing the excellent research facilities of 
the School of Physics. One of us (C.A.P.) is the 
holder of a University of Sydney Research Stu- 
dentship. 
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POSSIBILITY OF A 4-PION STRUCTURE IN THE T=J=1 PION-PION SCATTERING AMPLITUDE 
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Evidence has been accumulating of late on the 
existence of an enhanced or resonant energy re- 
gion in the T =J=1 pion-pion scattering ampli- 
tude. This evidence comes mainly from the 
analysis of the vector part of the nucleon form 
factors’ and from a fit of the low-energy pion- 
nucleon “small” phase shifts, in the combination 
fua-~fyo-** The qualitative agreement about the 
existence of such a resonant state has not been 
followed yet by a quantitative agreement about 
its location and width that, on the contrary, have 
been spanning a large range, depending on the 
authors and on the experiment analyzed. We wish 
to point out here that there could exist a structure 
effect, due to an enhanced 4-pion state, that 
could possibly account for the observed discrep- 
ancies and that could in principle be observed 
directly in production experiments. 

Let us analyze briefly how the many-pion con- 
figurations can be influenced by a resonance in 
the antisymmetric (both in isobaric spin and in 
configuration space) T=J=1 state. As we are 
mainly concerned with the possible contribution 
of such many-pion states to the nucleon form 
factors, we shall limit ourselves to the analysis 
of the many-pion “clouds” (JT =0 or 1).” The 
easiest way to describe the configuration of a 
“cloud” is by means of Young diagrams.® As 
the total wave function of the cloud has to be 
symmetric, we can limit ourselves to the anal- 
ysis of the configurations in isobaric space; an 
identical Young diagram will describe the state 
in configuration space. By using the notation in 
reference 7, the following results can be readily 
obtained: 

n=3: the maximum number of T =J=1 pion- 
pion bonds (“bi-pions”) is present in the config- 
uration (111), represented by Fig. 1(a). This 
configuration contains three bonds, and has T=0 
and J=1; it coincides with the enhanced 3-pion 
state discussed by Chew.® We shall call it, in 
what follows, the “three-pion.” 

n=4: the maximum number of “bi-pions” is 
present in the configuration (211), represented 
by Fig. 1(b). It contains three bonds, has T=1 
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FIG, 1. Young diagrams for some of the many-pion 
configurations. 


and J odd; it can be described as a “three-pion” 
plus a spare pion, which leads to an immediate 
understanding of its T=1. In its lowest angular 
momentum state it has the same quantum num- 
bers as the T=J=1 “bi-pion.” 

n=5: the maximum number of “bi-pions” is 
present in the configuration (221), represented 
by Fig. 1(c). It contains four bonds, and has 
T=1 andJ even. The configuration (311), which 
has T =0 and J odd, contains only three “bi-pions. 

It is clear from what precedes that should the 
“three-pion” state turn out to be strongly reso- 
nant or even bound, the 4-pion configuration (211) 
could play a determining role in the T =J=1 pion- 
pion scattering amplitude, where it would present 
itself as a second resonance or broad enhance- 
ment. Recent investigations on the scalar part 
of the nucleon form factors,* together with the 
results obtained by Abashian et al.,’° strongly 
suggest the T=0, J=1 configuration (“three-pion 
to be resonant at very low energy or, indeed, 
bound. It is also clear that the 5-pion (221) con- 
figuration, which can be described as a “three- 
pion” plus a “bi-pion,” cannot give trouble in the 
analysis of the vector part of the form factors, 
having an even J; while the (311) configuration, 
that could contribute to the scalar part of the form 
factors, is probably a much weaker configuration 

Having thus established that the existence of a 
“bi-pion” plus that of a resonating or bound 
“three-pion” makes the (211) T =J=1 “four-pion” 
an interesting object to investigate, we have tried 
to understand what information on the “four-pion 
could be inferred from the knowledge of the reso- 
nance (or binding) energies of the “bi-pion” and 
the “three-pion.” This is a quite standard prob- 
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lem in low-energy nuclear physics, where at 
least for small values of A one wants to be able 
to predict something about the ground state of 
nucleus A by using experimental information 
about the nucleus A -1 plus one’s knowledge of 
the forces between individual nucleons. It is 
quite a difficult and involved problem there, and 
it is even more difficult and involved in our case, 
of course, as we have no reliable “potential” 
which could describe the pion-pion interaction 
and therefore permit some sort of variational 
calculation for the 4-pion state. We have, how- 
ever, made an attempt in this direction by using 
the so-called projection theorem, in the form re- 
cently used by Iwao et al.,"" and applied by them 
to very light nuclei. We are fully aware of the 
fact that, since the “bi-pion” is only a resonance 
and not a bound state, the definition of an average 
kinetic energy for the infinite shell that we are 
filling with pions is, at best, a doubtful one. 
Anyway, we give our results here, with the 
hope that they could indicate at least qualitatively 
what is the influence of the “bi-pion” and the 
“three-pion” on the “four-pion.” It is clear that 
the (211) 4-pion configuration in its lowest angu- 
lar momentum state (J=1) can be most conven- 
iently described as some kind of hydrogen-like 
state of a pion in S wave with respect to a “three- 
pion.” If we now disregard the contribution (pos- 
itive or negative) to the binding energy of the sys- 
tem coming from this S-wave interaction, we are 
left with only two parameters, namely: the mass 
of the “bi-pion,” ty’, and that of the “three- 
pion,” to. The projection theorem applied to 
this simple model gives, for the mass of the 
“four -pion,” ty’™*: 
a - ao + to. (1) 

Let us make now the most optimistic assump- 
tion on the mass of the “three-pion,” namely, 
let us identify it with the “particle” observed 
in reference 10; this gives tg = 5u? and there- 
fore: 


t, = 22 y? 16 yp? 12 y? 9p 





t*= 38 yp? 30 p? 25 y? 19 yp? 


(2) 


It is clear from (2) that the 2-pion cut used in 
references 1-6 and approximated by a narrow 
Breit-Wigner formula or at least by a very en- 
hanced bump, could present on the contrary a 
well-defined structure (as the probable widths 


are smaller than the predicted separation of the 
bumps) and, what is more important, that the 
first of the two bumps (corresponding to the “bi- 
pion”) could also span a region corresponding to 
a much smaller fy than the values assumed at 
present. There is no real need for this enhance- 
ment to be really a very sharp resonance (see, 
for instance, reference 6); what matters here is 
that the “bi-pion” contribution could begin at 
lower energy and be followed by a “four-pion” 
contribution at higher energy. The actual de- 
termination of the parameters of an “equivalent” 
single resonance (up to now, attributed only to 
the “bi-pion”) could then be strongly sensitive 
to the method used for estimating its contribu- 
tion to the particular physical process studied, 
to the number of subtractions used in the repre- 
sentation of the physical amplitudes, and so on. 
This could perhaps offer one of the many possible 
justifications for the discrepancies found so far 
by the various authors in the determination of 

t 1-6 

"We note here that our neglect of the S-wave 
interaction can find some justification in the 
observational evidence for the absence of a (300) 
3-pion configuration which has T =1 and contains 
only symmetric 2-pion bonds.’° 

It is obvious that production experiments, anal- 

ogous to that by Abashian et al.’° with production 
of H® and at higher energy, could directly demon- 
strate the suggested two-bump behavior of the 
pion-pion T =J=1 scattering amplitude. Of 
course, this has nothing to do with the bump 
which is observed at present at a much lower 
energy, and which can be interpreted as the 
“three-pion” or, alternatively, as an S-wave 
pion-pion effect with T=0."* The second alterna- 
tive, that would leave only the form factor anal- 
ysis as an indication for the existence of the 
“three-pion,”* seems, however, in strong con- 
tradiction with the most recent experiments on 
production of pions by pions on protons,** where 
no indication of a T=0 (S-wave) state was found 
and where the evidence is for a strongly preferred 
T =1 production. Pion-production experiments, if 
analyzed by means of the extrapolation proce- 
dure," could indeed lead to a direct determina- 
tion of the pion-pion cross section; they contain, 
however, a mixture of T=0, 1, and 2 states (al- 
though charge ratios clearly indicate a preferred 
production of T =1 states) and a mixture of angu- 
lar momentum states, which makes the interpre- 
tation of the results not completely reliable. To 
this difficulty has to be added the well-known dif- 


73 














VoLUME 7, NUMBER 2 


PHYSICAL REVIEW LETTERS 





Jury 15, 1961 





ficulty of obtaining unambiguous answers from an 
extrapolation procedure. At best, this gives in- 
dications of the general features of the process 
analyzed; to get more refined details and struc- 
ture indications from the experimental data would 
require much more complete experiments than 
those available at present. 

Useful conversations with Professor M. Cini 
are gratefully acknowledged. One of us (V.D.A.) 
wants also to thank warmly the Istituto di Fisica 
di Roma for the hospitality extended to him. 





‘w. Frazer and J. Fulco, Phys. Rev. Letters 2, 
365 (1959). 

*wW. Frazer and J. Fulco, Phys. Rev. 117, 1609 
(1960). 

4S, Bergia, A. Stanghellini, S. Fubini, and C. Villi, 
Phys. Rev. Letters 6, 367 (1961). 

‘F, J. Bowcock, W. N. Cottingham, and D. Lurié, 
Nuovo cimento 16, 918 (1960); Nuovo cimento 19, 142 
(1961); Phys. Rev. Letters 5, 386 (1960). 









5s, Frautschi, Phys. Rev. Letters 5, 159 (1960). 
6J, Hamilton, P. Menotti, T. D. Spearman, and 
W. S. Woolcock, Nuovo cimento 20, 519 (1961). 

TA, Pais, Ann. Phys. 9, 548 (1960). 

SYoung diagrams were first applied to the study of 
the 3-pion configuration by D. Amati and B. Vitale, 
Nuovo cimento 2, 719 (1955); a very extensive appli- 
cation to the general case of n-pion “clouds” has been 
given by Pais, reference 7. 

9G. Chew, Phys. Rev. Letters 4, 142 (1960). 

104, Abashian, N. E. Booth, and K. M. Crowe, Phys, 
Rev. Letters 5, 258 (1960). 

11g, Iwao, F. B. Wang, and F. Morin, Atomic Energy 
Commission Report NYO-9536 (to be published). 

27, N. Truong, Phys. Rev. Letters 6, 308 (1961). 
133, A, Anderson, V. X. Bang, G. Burke, D. D, 
Carmony, and N. Schmitz, Phys. Rev. Letters 6, 365 
(1961); A. R. Erwin, R. March, W. D. Walker, and 
E. West, Phys. Rev. Letters 6, 628 (1961); D. Stone- 
hill, C. Baltay, H. Courant, W. Fickinger, E. C. Fow- 
ler, H. Kraybill, J. Sandweiss, J. Sanford, and H. Taft, 
Phys. Rev. Letters 6, 624 (1961). 

“4G. Chew and F. Low, Phys. Rev. 113, 1640 (1959). 





